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Abstract. We suggest a novel proposal to express decoherence in open quantum 
systems by jointly employing spectral and stochastic methods. This proposal, which 
basically perturbs the unitary evolution operator in a random fashion, allows us to 
embrace both markovian and nonmarkovian situations with little extra effort. We argue 
that it can be very suitable to deal with models where an approximation neglecting 
some degrees of freedom is undertaken. Mathematical simplicity is also obtained both 
to solve some master equations and to arrive at experimentally measured decoherence 
functions. 
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1. Introduction 

Decoherence in an open quantum system is the loss of quantum coherence due to 
the entanglement with its environment as a consequence of their mutual interaction 
It precludes thus the formation of some superpositions of states, creating in the 
system environment-induced superselection rules p], [3|. The rigorous mathematical 
expression of this phenomenon runs parallel to this scheme, namely, one should first 
consider the joint evolution of both the system and the environment and then take 
the partial trace over the environmental degrees of freedom || ||. In general, the 
resulting evolution equation for the density operator of the system ps is non-markovian 
and extremely difficult to handle. This is the reason why suitable approximations have 
been traditionally pursued. In particular, adequate methods to derive an approximate 
markovian equation have been found |5], [| |7|. On the other hand, more axiomatic 
results have also been developed which pose some general physical conditions for the 
evolution of the system density operator || ^] and arrive at the general form for the 
evolution equation of ps- Both approaches show pros and cons: the axiomatic program 
starts by assuming the markovianity of the evolution as one of the physical hypotheses, 
thus making rather difficult the generalization to non-markovian situations, which, due 
to experimental progress, are becoming interesting enough to develop theoretical tools. 
On the other hand, the constructive approach depends very sensitively on the particular 
model of interaction between the system and its environment, thus rendering it rather 
inappropiate for generic discussion beyond precise details. 

Here we suggest another proposal which embraces both markovian and non-markovian 
situations, like the constructive approach, but it does not depend, in general terms, 
on the specific features of the system-environment interaction. The physical picture 
which motivates this scheme models the effect of the environment upon the system by 
a random perturbation of the evolution operator. Since the origin of the stochasticity 
is rooted upon the environment and since we do not have control over its degrees of 
freedom, we then take the stochastic expectation value. This is of course not new, 



and indeed it has been used in classical physics too [10, 11]. The novel proposal stems 



from the combination of this old idea with the spectral decomposition theorem for the 



evolution operator Jl2], [L3|] . This combination will turn out to be quite useful to express 
decoherence in certain models. 

The paper is organized as follows. In section ^| the general scheme of this proposal is 
presented, then in section |3| the case of global random perturbations of the Hamiltonian 
spectrum is studied to arrive at a Lindblad-type evolution equation both when the 
Lindblad operators are selfadjoint (subsection |3.1| ) and nonself adjoint (subsection |3.2| ). 
In the following section (|J) a detailed example of the conjuction of the spectral 
decomposition and stochastic methods is worked out. A brief discussion and some 
conclusions are then presented in section |5|. Finally we close the paper with a 



computational appendix ( |Appendix Aj) and a mathematical proof ( Appendix B|) . 
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2. General scheme 

To show how stochastic methods and the spectral decomposition theorem are combined, 
consider a quantum system with hamiltonian H, which for simplicity's sake will be 
discrete and nondegenerate. The evolution equation for this system will then be given by 
the evolution operator U (t) = exp(-itH), which by means of the spectral decomposition 
of H = e kPk can be expressed as (ft = 1 as usual) 

U(t) = Y,^ iekt Pk (2.1) 

k 

The operator U(t) will then be randomly perturbed as the following relation shows 
U ® = E e ~' lektPk ^ U «® = E e ' iXk(t)p k (2-2) 

k k 

where Xk(t) denotes a real- valued stochastic process for each k. The system density 
operator ps(t) at time t will then be given by the expectation value ps(t) = 
^[U s t(t)p(0)Ug t (t)]. The density operator, in the Hamiltonian eigenvector basis, is then 
given by 

p n , m (t) = E[ e -^W-x«( t »]p him (0) (2.3) 

where it has been supposed that the system and the environment are not initially 
correlated, thus p(0) is not random. Both the original unitary evolution and the 
decoherence function are contained in the expectation value Wye~ liyXniyt ^ Xmiyt ^} and should 
then be properly disentangled. To do this note that the "stochastic promotion" 
expressed by ( |2.2|) can be very intuitively understood if we express the argument of 
the exponential in the original evolution operator as an integral and the stochastic 
process Xk(t) is written as an Ito process fH 

e k t= / e k ds -> Xk(t) = / a k {s)ds + / b k (s)dB k (s) (2.4) 
Jo Jo Jo 

where a k (t) and b k (t) are real- valued stochastic processes for each k and B k (t) are 
standard real Brownian motions for each k. 

Under full generality the previous expectation value cannot be computed, but with 
some general physical assumptions some useful results may rapidly arise. As a very 
elementary example consider for instance the white-noise perturbation, i.e. 

e k t -> e k t + / e k £(s)ds = e k t + e k B k (t) (2.5) 
Jo 

where £(£) denotes white noise with correlation properties E[£(£)] = 0, E[£(i)£(s)] = 
j8(t — s). Note that all eigenvalues are perturbed by the same stochastic noise £(t). 
Under this assumption and after computing the expectation value 



E[g- i (efc-e fc /)/ t C(s)dsj __ e ~^(e fc -e fe ,) 2 
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the elements of ps{t) in the Hamiltonian eigenvector basis are given by 

PnAt) = P^)^"^ (2-6) 

where p^\t) denotes the unitarily evolved density operator. The corresponding 
evolution equation for ps is straightforwardly shown to be the markovian phase- 
destroying master equation (see [JTH] for an alternative deduction) 



^ = -i[H,p s (t)} - 1[H, [H,p s (t)\] (2.7) 

Thus the effect of the environment can be understood as a global random energy 
"kick" upon the system. This method also provides a very fast way to solve equation 
( P-7| ) provided one knows the solution to the unitary evolution (cf. |16| for an explicit 
example). 

3. Global perturbations 

3.1. Self adjoint Lindblads 

The global random perturbations like in the example in the preceding section can 
also be expressed without resorting to the spectral decomposition theorem. A partial 



treatment was provided in [jig] , where the most general Lindbladian evolution with 
selfadjoint Lindblad operators was obtained. Here we offer a complementary proof of this 
result which uses different techniques and thus paves the way for later generalizations 
(nonself adjoint Lindblad operators). Let H be the Hamiltonian of the system. Introduce 
now a global random perturbation by adding the operator r](t)V to H, where rj(t) 
denotes an arbitrary real-valued stochastic process and V an arbitrary selfadjoint 
operator. The evolution operator will then read 

U st {t) = exp[-itH-i n {t)V] (3.1) 

The main trouble to compute p(t) = E[U st (t)p(0)Ul t (t)] stems from the 



commutativity properties of H and V. If [H, V) — it can be readily shown [|T^ 
that the evolution equation for ps is 

^ = -i[H,p s (t)} - ~(vW W, /*(*)]] (3-2) 

where (r]) t denotes the quadratic variation of r](t) [|17|]. When [H, V] ^ 0, the 
computation cannot be performed in the same way. An option to overcome this difficulty 



was presented in ||16|| , where the stochastic perturbation was performed in the Heisenberg 
picture; here we resort to Feynman's operational calculus |18|, [19], |20| . In particular, 
the following three Feynman's heuristic rules to deal with functions of noncommuting 
operators will be thoroughly used: 



(i) Attach time indices to the operators to specify the order of operators in products. 



Expressing Decoherence with Spectral and Stochastic Methods 



5 



(ii) With time indices attached, form functions of these operators by treating them as 
though they were commuting. 

(iii) Finally, "disentangle" the resulting expressions; that is, restore the conventional 
ordering of the operators. 



A mathematical formalization of each of these rules can be found in [[B], [ZIJ. Here 
we will only show how they can be used to prove our result. Consider then the evolution 
operator ( [3.1|) ; attach time indices both to H and V so that one can write 



U st (t) = exp f-ijf H{s)ds-\j V(s)dfi(s)^J 



(3.3a) 



where dfi(t) denotes the Lebesgue-Stieltjes measure rj(t)dt. Now treat both factors as 
though they were commuting: 



(3.3b) 



U s t(t) = exp i J H(s)ds ) exp ( — i / V(s)dfi(s] 
The density operator then reads 



p{t) = e-i/o^W^E 



exp -i / V(s)dfj,(s) p(0) exp i / V(s)dfi(s 



3^o^ (s)ds (3.3c) 



Now the expectation value can be computed using (see also appendix |Appendix A| ) 
e -iA Be iA = e -i[A,-) B with A = A t > w here e~^B = £~ ^[A, [A, [A, [A, B]]]]: 



E 



exp ( — i / V(s)dfj,(s) I p(0) exp ( i / V(s)dfi(s 



exp - 



V\s)d{ri) t 



p(0)(3.3d) 



And finally, restore the conventional ordering of the operators: 
Ps(t) = exp (-i[H,.]-~(ri) t [V 2 r }\ p(0) 
which drives us to the master equation 



(3.3e) 



(3.3/) 



We have proved the main result found in JTB| with other means: Any Lindbladian 
master equation, either markovian or non-markovian, with selfadjoint Lindblad 
operators, can be understood as a random unitary evolution. The generalization to many 
Lindblad operators can be readily accomplished by using several uncorrelated stochastic 
processes r)i(t) i — 1, . . . , n. Note that the markovian case is obtained only when (i]) t is 
linear in time, i.e. when rj(t) = const. +7 1//2 £> t which amounts to perturbing the original 
Hamiltonian by a white noise like in (|2.5|). 



Complementarily one can also use these techniques to study the evolution of a 
system given by a Hamiltonian with one or several stochastic parameters in it. Take 
for instance the system with Hamiltonian H = H\ + X(t)H 2 and let the parameter 
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X(t) be uncontrollably random. Following the previous prescriptions the corresponding 
evolution operator will read 



U st {t) = Texp 



itHi - \H 2 //(*) + / \{s)dB{s 



(3.4) 



where T denotes the time-ordering operator and by definition J* X(s)ds = (j,(t) + 
Jq X(s)dB(s). The master equation for ps will then be 

^ps(t) = -i[H llPs (t)]-im[H2,Ps(t)}-~Mt)+ f ~X(s)dB(s)) t [H 2l [# 2 ,p s (t)]](3.5) 



dt 



2dt 



This scheme has already been used pi] , pT6|| to study the effects of laser intensity 
and phase fluctuations in ion traps. In more general terms, these techniques appear 
very suitable to deal with approximate models in which the approximation renders the 
physical realistic situation mathematically manageable to perform computations, but on 
the other hand it ineludibly carries some decohering effect in it . An explicit example was 
worked out in ||16|| , where the Jaynes-Cummings model (JCM hereafter) was revisited 
within this approach. The JCM is a crude, though very appropiate, approximation 
to deal with the interaction between an atom and the electromagnetic field in certain 
circumstances (cf. |22], |23| ) which reduces both systems to a two-level system and a 
harmonic oscillator, respectively, interacting via a Hamiltonian Hj = A(cr_a^ + o~ + a), 
where A is the coupling constant, a± are the raising/lowering atomic operators and 
a(a^) are the annihilation (creation) field mode operators. It is then clear that the 
rest of field modes as well as other atomic degrees of freedom are completely neglected, 
thus providing an ideal situation to quantitatively describe some physical phenomena 
like Rabi oscillations, periodic collapses and revivals or squeezing [p4 |. Decoherence is 



absent in this model, however experimental results |25], |26[ show the opposite, though 
extremely difficult to be attributed to the environment. This tiny decoherence effect 
may be incorporated into the model from the beginning by letting the coupling constant 
A have a random character, this randomness coming from the effect of the neglected 
degrees of freedom. The theoretical predictions are in good agreement with experimental 
results ||16|| . Note that within this approach decoherence should be understood as an 
intrinsic phenomenon stemming out from the adopted approximation and thus being 
an inherent property in those models which neglects degrees of freedom of the system. 
Notice how despite providing the same analytical behaviour, this approach conceptually 
differs from others in which a quantum bosonic reservoir is explicitly taken into account 
27L |28 . 



3.2. Nonself adjoint Lindblads 

A natural question arises concerning whether it is possible or not ot arrive at a 
Lindbladian master equation with Lindblad operators not necessarily selfadjoint like 
in the previous examples. And the answer is positive but with important remarks. 
To be concrete let us consider a spin-1/2 interacting with a random magnetic field 
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B(£) = B x (t)~x + B y (t)y + 5 z (£)x. The Hamiltonian describing the evolution of the spin 
will then be 

H(t) = lj o- z + B(t) • er (3.6) 

Complementarily this can be viewed as a quantum bit (left term of the rhs of ( |3.6D ) 
undertaking a logical operation in a noisy environment (right term of the rhs). Without 
loss of generality we can focus on those cases in which B(£) • z = 0, since otherwise we 
can use the tools developed in the preceding section. The evolution operator will then 
be 



Ust(t) = Texp 



-iu t — i (B x (s)a x + B y (s)a y ) ds 



(3.7a) 



Since we will make use of Feynman's operational calculus as in previous paragraphs, 
we can concentrate on the operator 



Ust(t) = Texp 



i/ (B x (s)a x (s) + B y (s)<7y(s)) ds 



(3.76) 



The magnetic field components will be expressed through Ito integrals as 



t rt 

B x (s)ds = / b x (s)dB x (s) (3.8a) 
o Jo 

B y (s)ds = / b y (s)dB y {s) (3.8b) 
o Jo 

where the b k (tys are arbitrary deterministic real- valued functions and Bk(t) denotes 
nonstandard real Brownian motions (cf. below). The careful reader may argue at first 
sight that this is not the more general form for a real stochastic process, and it is not, 
though it contains all physically relevant cases which interest us in this analysis. To 
prove that, we begin by demanding that the magnetic field components be martingales 



14, 17], which is general enough from a physical standpoint. Then making use of the 



martingale representation theorem, they can be uniquely written as 



B k (s)ds = E 



o 



B k (s)ds 



g k (s)dB k (s) k = x,y (3.9) 



where g k (t) is a real stochastic process depending on the particular choice of B k (t). 
Without loss of generality we take E[Sjt] = since otherwise this factor can be included 
in the dropped-out exponential after application of Feynman's rules and finally recovered 
back in the final step. We are interested in those effects coming out from the stochasticity 
of the Hamiltonian and which are at the end detected after calculating the stochastic 
average. This is the reason why we can also consider the processes (?&(£) 's as being 
deterministic. 

The correlations properties of the magnetic field components will be generically 
expressed through the correlation properties of the Brownian motions 
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E[B 2 x (t)} = lx t 
E[B 2 y (t)] = ly t 
E[B x (t)B y (t)] = 7xy t 



(3.10a) 
(3.106) 
(3.10c) 



In order to expouse clearly the role played by the correlation properties between 
both componentes, let us concentrate upon white-noise process, i.e. 6 x (i) = b y (t) = 1 
(equivalently B x (t) = B y {t) = £(t)). Now we decompose a x and a y in terms of the 
raising and lowering operators to rewrite (|3.7fl) as 



Ust(t) = Texp 



— i 



[B(s)a+(s) + B*(s)a^(s))ds 



(3.11a) 



where B(t) is a complex stochastic process with %KB(t) = B x (t) and 3B(t) = —B y (t). 
Alternatively equation (|3.11o| ) can be rewritten as 



Ust(t) = Texp 



-i / <r+(s)dB(s) - i / c-(s)dB*(s) 



(3.116) 



where B(t) is a complex stochastic process with correlation properties given by 



E[B 2 (t)] = ( 7x - 7y - i2 7xy ) t (3.12a) 
E[B(t)B*(t)] = ( lx + ly )t (3.126) 

Now using again e~ lA Be lA = e~^ A '^B with A — we compute the stochastic 



expectation value p(t) = E U st {t)p{Q)U i 



'st 



pit) = TE 



[p(0)](3.13) 



exp f-ij [<j + (s) r ]dB(s)-ij [<7_(s),-]d6*(s) 

After computing this expectation value (see appendix |Appcndix A ) and restoring 
the conventional time ordering of operators, the final result is (up to Hamiltonian factors) 



where 



p{t) = exp [t (£ xy + £+ + £_+ V xy )} [p(0)] 



(3.14) 



XV 



V 



XV 



7x-7y 



{Wx, [0"x,']] 



(7x + 7y) (h-,^] + h,^/] 

2 7 xy (0"x • Oy + 0"y ' Ox) 



J = +, 



(3.15a) 

(3.156) 
(3.15c) 



The first two generators are Lindblad generators, thus denoting a completely 
positive markovian evolution || |J whereas the last one is not. Since the evolution 
is clearly markovian, this means that complete positivity is not ensured during the 
evolution. It is remarkable how once more the presence of correlations spoils the 
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complete positivity |29], |3(| . In fact this is a clear example of how correlations within 



the environment, even in the case of no correlations between the system and this 
environment, prevents the evolution from being completely positive (see Appendix B for 



a mathematical proof of this fact). Notice however that when the random magnetic field 
is isotropic (7* = 7 y ) and the components are uncorrelated (7 xy = 0), a fully Lindblad 
evolution is recovered with dissipative generator given by Q3.156| ), i.e. a Lindblad 
evolution with nonself adjoint Lindblad operators. Nevertheless this generator shows 
two opposite parts which within this approach cannot be separated. Physically this is 
rooted upon the isotropy of the magnetic field which, as expected, drives the system 
with equal probability both to the ground and to the excited state. 
The generalization to nonmarkovian evolution can be performed with little extra effort. 
For simplicity we keep the isotropy condition, thus 6 x (i) = b y (t) = b(t). The random 
evolution operator will then be 

U st (t) = Texp ^-i b(s)a+(s)dB(s) - i jf b(s)o-_(s)dB*(s)^J (3.16) 

where the complex stochastic process B(t) have the same correlation properties as above. 
The calculation proceeds along similar lines with result 

pit) = exp [A(t) (£ xy + £+ + £_ + V xy )] [p(0)) (3.17) 

where the generators are the same as before and A(t) = j^b 2 {s)ds. This is clearly 
nonmarkovian and reduces to the markovian case only when b(t) = const. Regretfully 
the previous calculations show that the dissipative generator appears always in opposite 
pairs. The conditions under which a single dissipative Lindblad generator is obtained 
are being currently under study (see also [RTI] for an alternative approach). 



4. Local perturbations 

But the versatility of the spectral decomposition theorem appears when the perturbation 
does not affect the whole system globally, as in the previous examples. By this we mean 
that in the random perturbation (|2.4j ) the Brownian motions Bj.(t) do not show singular 
correlation properties!], but rather general ones. Take for instance the case in which 

e k t -> e k t + >ye k B k (t) (4.1) 

with correlation properties E[B k (t)B k >(t)} = 5 kk >t. Under this assumption (which 
physically can be understood as the independent or incoherent effect of the environment 
upon each system energy level), the density operator components are given by 



Pnm{t) = pJ2(t) n = m (4.2a) 
= P.|lWe"^ (4+4l) n + m (4.25) 

I Note that in the previous cases -global perturbations- the processes Bk(t)'s satisfied E{Bk(t)Bk> (t)] = 
t for all pairs k, k', i.e. all Bk(t)'s were esentially copies of the same standard one B(t) in (2.5). 
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which show a faster decoherence than ( f2.6| ). By choosing another correlation properties 
among the Bfc(t)'s one can arrive at more general decoherence functions for the decay 
of off-diagonal terms of the density operator. As a prominent example let us consider 
two quantum particles interacting through a central potential V(r 12 ). The Hamiltonian 
will then be given by 

H = H° 1 + H° + H l2 = -g^Ai - ^-A 2 + AV(r 12 ) (4.3) 

where A is a coupling constant and r 12 = |r x — r 2 |. In terms of the center-of-mass 
Hcm = R- an d the relative position r 12 = r coordinates, the Hamiltonian is given by 

with M and fi the total and reduced masses respectively. The total energy is written as 
E T = Ecu + E rel (4.5) 
so that the corresponding evolution operator will then be given by 

U(t) = [ e-' ltET dP(E T ) (4.6) 

Now a weak environment-system coupling can be expressed through a random 
perturbation of the total energy Et by a stochastic amount generically expressed like 
J e(s; E T )dB(s; E T ) where e(s; E T ) = e(s; E C m) only depends on the center of mass and 
where the standard real Brownian motions show correlation properties only dependent 
on the center of mass variables as well. This is the main physical hypothesis and it is 
expressed as follows 

U st (t) = j e -itE T -i^e{s;E T )dB{s;E T ) dp ^ ^ 

= / e -itE T -iJ t e(s;EcM)dB(s;E T ) ( j t p^ T ^ (4-7&) 

with E[B(t; Et)B(J:; E' t )] = g(t; Et, E' t ). Note that the correlation function satisfies 



g{t;E T ,E T ) = 1 (4.8a) 

g(t;E T ,E' T ) < 1 (4.85) 

g(t, E T , E' T ) = g(t; E CM , E' CM ) (4.8c) 

where the last equation reflects the assumed physical hypothesis. Then if we denote the 
stochastic process 

rj(t; E T , E' t ) = [ e{s; E T )dB{s; E T ) - [ e(s; E' T )dB{s; E' T ) 

Jo Jo 
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the components of the density operator in the basis [Ecu] E ie i) will then be given by 

(E CM] E reX \p(t)\E' CM ; E' Iel ) = P (0) (^cm, E iel ; E' GM , i^-fW-;^))* (4.9) 

Now the quadratic variation of n can be calculated in terms of the correlation 
function g above: 

M-; E T , E' T )) t = [ [e 2 {s; E T ) + e 2 {s; E' T ) - 2g{s; E T , E' T )e(s; E T )e{s; E' T )] ds = 
Jo 

= / [e 2 (s;E CM ) + e 2 (s;E' CM ) - 2g(s; Ecu, E' cu )e(s; E C M)e(s; E' cu )] ds 
Jo 

= (fj{-,E CM ,E' CM )) t (4.10) 

Note that (?7(-; Et, -Er))t = 0, thus the trace is conserved. Also notice 
that components showing quantum coherence for different center-of-mass modes 
asymptotically vanish, since 

^(ry(-; E T , E' T )) t > (e(s; E CM ) - e(s; E' CM )) 2 > (4.11) 

On the other hand, full quantum coherence is kept for the relative coordinate, thus 
internal dynamics is not affected. A particular example of this scheme has been applied 
to arrive at the experimental expressions obtained in the study of Rabi oscillations in 
an ion trap Jl6| . 

This same technique can also be applied in those cases in which the randomness appears 
only in one part of the total Hamiltonian. Let for instance H = H + Hi be the total 
Hamiltonian of a quantum system and suppose that the randomness affects only the 
interaction part of H, i.e. Hj. One has two complementary alternatives. On one 
hand, one may change to the ifo-picture and then apply the preceding techniques to 
the evolution operator Ui nt (t) = Texp ^— i Hi(s)dsj where H\{t) denotes the H\ 

Hamiltonian in H representation, i.e. Hi(t) = e °Hie~ ltH °. On the other hand, one 
may also resort to Feynman's operational calculus and "break" the exponential into two 
parts: a first deterministic exponential involving only Ho and a second one involving the 
random interaction Hamiltonian Hj, which can now be expressed through the spectral 
decompostion theorem and then stochasticly promoted like above. Both approaches are 
equivalent. 



5. Discussion and conclusions 

The use of stochastic methods in Hilbert space is of course not new (cf. |5|, Q 
and references therein). Here they have been used in conjuction with the spectral 
decomposition to accomodate a great variety of decoherence functions which may appear 
in realistic situations and which otherwise require more or less complex calculations. It 
is important to remark that the whole scheme presented above preserves the unitarity 
of the evolution operator, i.e. U s t(t) is unitary with probability 1. This is in contrast to 
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the radically different usage of these methods in stochastic dynamical reduction theories 
(cf. |33| and references therein for diverse examples within this philosophy), where the 
norm of the initial state vector is not preserved (reduction) and thus one is obliged to 
introduce non-linear terms to keep the probability interpretation. A clear-cut example 
of the difference between the decoherence and the reduction behaviours appears in |34[ . 
Possible connections between these two distinct pictures (decoherence and reduction) 
are currently under study We argue that these possible connections, if any, appear 
in those cases where an analytical continuation in the complex plane of the spectrum 
of the Hamiltonian is used ||35|| , which requires extending the quantum mathematical 
formalism of Hilbert spaces to that of rigged Hilbert spaces. We further conjecture 
that obtaining a Lindbladian evolution equation with just one nonselfadjoint Lindblad 
operator is also related to the existence of complex Hamiltonian eigenvalues and thus 
to micro- irreversibility [|35| . 

In author's opinion, a remarkable advantage of the approach presented above (apart 
from some mathematical benefits in solving some master equations |15|, pl|) stems from 
the possibility of dealing with approximate (though necessary) descriptions such as the 
JCM, where the approximation itself carries an inherent decohering effect as a result 
of neglecting some degrees of freedom of the system. This effect has been expressed 
through the randomization of some internal parameters of the system Hamiltonian and 
contrasted with some experimental data |16|]. The good agreement achieved invites us 
to look for new realistic situations in which we can further check the utility of this 
approach. A second advantage obtained from the involved mathematical formalism is 
the easy and fast generalization to nonmarkovian situations. However the approach 
is of phenomenological nature, consequently if one is interested in specific details 
depending e.g. on the concrete interaction between the system and the environment, it 
is compulsory to resort to the well-known "tracing-out" methods 

To sum up, in this work we have shown how the conjuction of stochastic methods 
and the spectral decomposition theorem can provide the analytical expression for the 
experimental behaviour of certain decohering quantum systems in a rather simple form. 
The main idea consists of identifying the pervasive effect of the surrounding environment 
upon a quantum system as a random perturbation of its evolution operator. Diverse 
decoherence functions can be obtained by adequately choosing the correlation properties 
of the random parameters. As a matter of fact it can be proved that when there 
exists correlation between these different parameters, then the system evolution is not 
completely positive, providing a further physical condition by which complete positivity 
is broken. It has been shown, as a noticeable example, how in a compound system the 
center of mass may suffer from decohering effects whereas the internal dynamics, i.e. the 
internal modes of the system are completely unaffected. The utility of this approach to 
deal with some approximate descriptions of quantum systems has been briefly discussed. 
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Appendix A. Computation of E [exp (— i2l(£))] 

Here we include the computation of the stochastic average of the exponential of a 
random superoperator 2l(t) = f* /3(s)&(s)dB(s), where /3(s) is a deterministic function, 
&(s) is an arbitrary superoperator and B(s) is a complex or real Brownian motion. 
By linearity of the expectation value, the computation reduces to find the nth order 
moment of To do this, notice that it can be expressed through the Ito equation 

dSL(t) = P(t)&(t)dB(t) and then applying Ito's formula to the function ^(2t(t)) = 2l n (t). 
The nth order moment 9Jt n (£) = E[2l n (t)] then satisfies the recursive differential equation 

dWl n (t) = n{n ~ l) VR n „ 2 (t)p 2 (t)& 2 (t)d{B) t n > 2 (A.l) 
Recalling that Wl = 1 and DJti = one finally arrives at 



m n+ i(t)=0 (A.3) 

where &(t) = f* (3 2 (s)& 2 (s)d(B) s . Substituting these relations in the series development 
of the exponential one immediately arrives at 

E[exp(-i2l(t))] = expf-^ (A.4) 



Appendix B. Initially correlated environments might preclude completely 
positive system evolutions 

Here we prove how the presence of quantum correlations inside the environment might 
prevent the system evolution from being completely positive, even when the system is 
not initially correlated with the environment. Let the whole system+environment be 
represented by a Hilbert space <8> £h <8> #2, where corresponds to the system and 
<8> #2 to the environment. Let the respective dimensions be N$, Aq and N 2 and 
the generators of the group SU(Nk) be denoted by a^f \ with j — 1, . . . , N k . The most 
general joint density operator can then be written as 



PS12 



N S N ± N 2 



hu + ctio-f ) <g> Ii <g> I 2 + Pjh <E> ® I 2 + 7fcls <E> Ii <8> o-f ) 



+ ^(jf } <g> <g> I 2 + e ik <jf ) <g> Ii (8) ai 2) + r^Is ® cr) 1} ® a[ 2] + 

+ u ijk af ] <g> ® <r< 2) ] (B. 

Then denoting p fe (£) = Tri m (u S i2(t)psi2(0)l4 12 (t) j with / ^ k ^ m, it is 
elementary to prove that psi2 can be rewritten as 
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where 



PS12 = PS ® Pi ® P2 + ) ® ® I 2 + 4^i S ] ® II ® + 

+ ^ fc I s ® of ® af + 3 ® of ® af (B.2) 



^"ivyvyv^^-^ AsaW (% " aaj) (2 - 3a) 



1 

NsN x N 2 V '" J ™ '■" ' lJK ~ N s N~iN 2 



= lulvnvT (%• - &7i) 4fe = IyHyHyT ( v W - OiPjl/k) (2-36) 



From this and after well-known arrangements [29], it is easy to arrive at 



Ps{t) = E M ^,a'6'(*)Ps(0)M a t 6)a , 6 , + 

ab a' 6' 

+ <^Tr 12 (U(t)(jf } ® of ® I 2 Ut(t)) + e; k Tr 12 (u(t)af > ® I x ® of ^(t; 



+ ^.Tr 12 (u(t)I s <g> of ® o£V(t)) + ^ k Tr 12 (u(t)of ) ® a} 1 ) <g) erf ljt(t) 

(2.4) 



From ( |2.4j ) it is clear that even when the system and the environment are not 
initially correlated ■ = e' ik = v'^ k = 0), the initial correlations inside the environment 
7)'j k 7^ might introduce a non-Kraus term in the evolved reduced density matrix of the 
system, thus precluding complete positivity. An example has been provided in the text 
using different methods. 

However note that when the coefficients satisfy 

«i = c\i (3j = fa 7 fc = 7fc (2.5 a) 

8ij = ai/3j e ik = aa k rjjk = fjjk (2.56) 
Vijk = aif} jk (2.5 c) 

for arbitrary "barred" quantities, the joint density operator can then be factorized as 
Psi2 = Ps <8> P12 and the reduced dynamics is then completely positive as expected. 
Another sufficient condition to get a completely positive evolution even in the case of 
initial correlation arises when the dynamics is not entangling |[36||, i.e. when the evolution 



operator can be written as Usi2(t) = Us(t) <8> U\{t) ® U2(t) (see for details). 
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